Correlation measurement of squeezed light 
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We study the implementation of a correlation measurement technique for the characterization of 
squeezed light which is nearly free of electronic noise. With two different sources of squeezed light, we 
show that the sign of the covariance coefficient, revealed from the time resolved correlation data, is 
witnessing the presence of squeezing in the system. Furthermore, we estimate the degree of squeezing 
using the correlation method and compare it to the standard homodyne measurement scheme. We 
show that the role of electronic detector noise is minimized using the correlation approach as opposed 
to homodyning where it often becomes a crucial issue. 
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I. INTRODUCTION 

Since thepioneering experiments of Hanbury-Brown 
and Twiss T] (HBT) on the realization of the stellar in- 
tensity interferometer, the correlation measurement tech- 
nique has proven its great potential in modern quantum 
optics. A direct observation of the second order intensity 
correlation function reveals many nontrivial properties 
of the photon statistics [1]. For example, the observa- 
tion of the strong photon bunching is common in exper- 
iments where correlated photon pairs are generated in 
various nonlinear optical processes. In practice, such a 
technique is involved in many fundamental and practical 
tasks such as tests of the foundation of quantum mechan- 
ics [1| , quantum cryptography [J , quantum metrology Q 
and many others. 

In contrast, in experiments with squeezed light, nor- 
mally associated with the detection of a continuous de- 
gree of freedom, a different measurement technique is 
traditionally implemented. The conventional method ap- 
plied in this case is balanced homodyne detection (HD) 
where the state under investigation (often a squeezed 
state) is mixed on a symmetric beamsplitter with a local 
oscillator (LO) and the difference photocurrent of two 
analogue photodetectors is recorded [1] (see also Fig.l). 
By scanning the phase of the LO the statistics of the 
field's quadratures are probed which subsequently allows 
for the full tomographic characterization of the state ■ 

Recently, there has been a trend to combine the mea- 
surement techniques of the discrete and continuous vari- 
able domains. In ref. HD techniques were applied to 
the characterization of a single photon Fock state as an 
alternative to the traditional correlation method. More- 
over the characterization of a squeezed state through 
measurements of the second-order correlation function - 
commonly used in the discrete variable quantum optics 
- was theoretically proposed in and experimentally 
demonstrated using single photon detectors [ll| and ho- 



modyne detectors [T2| . 

In the appendix of ref. p^ . the correlation measure- 
ment technique of squeezed light was analyzed and it was 
found that it provides an electronic-noise free approach 
for estimating the degree of squeezing. This electronic- 
noise free estimation method was experimentally imple- 
mented in ref. ^15.] for a bright squeezed light beam gen- 
erated in a fiber exploiting the optical Kerr effect. In 
the present paper we extend the experimental analysis 
to include vacuum squeezed light produced by an op- 
tical parametric oscillator. The utilization of squeezed 
vacuum yields a more thorough analysis of the effect of 
electronic noise on the measurement of squeezing. 



II. THEORY 

Let us consider a HD setup as shown in Fig.l: The 
signal field under investigation is mixed with a LO on 
a 50/50 beam splitter. The two outputs are measured 
with analogue detectors and the difference of the pho- 
tocurrents of the two detectors is produced. The theo- 
retical treatment is straightforward and follows the stan- 
dard formalism [31. The field operator can be written 
as a — a -\- 5a, where a represents the classical bright 
component and 5a is an operator with zero mean value, 
describing the quantum uncertainty of the field ampli- 
tude. We consider the input state to be a squeezed vac- 
uum state. Assuming unity detector efficiency we calcu- 
late the photocurrents of the two detectors in each port 
of the beamsplitter: 

ii,2 = l/2aio + ctLo{5XLo ± 5X^) -I- 5ieii,2 (1) 

where a^o and (j) is the mean field amplitude and the 
phase of the LO, respectively, 5X^ = cos 4)5 X + sin^JF 
is the fluctuation of the rotated quadrature of the in- 
put beam, SX^o is the amplitude quadrature of the LO 
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and Sieii.2 is a stochastic variable accounting for the elec- 
tronic noise (EN). To complete the homodyne setup the 
difference photocurrent is recorded and its variance is 
given by 

- = Aalo{5Xl) + {5^l,) + {Sr^^). (2) 

which is clearly seen to be measuring the quadrature fluc- 
tuations of the signal, (i5X|). However, it is also clear 
that the contribution from the EN always affects the 
measurement: The ENs are additive and their presence 
might dominate the contribution from the signal thus po- 
tentially masking the quantum features. The EN relative 
to the signal noise can be made very small partly through 
a clever design of the amplifying circuit of the detectors 
and partly by using a very strong LO which "amplifies" 
the signal contribution but leaving the EN unchanged. 
This is however not always possible and the performance 
of the HD will be ultimately limited by the EN. 

For the estimation of squeezing, a careful calibration 
of the shot noise level (SNL) must be carried out. The 
SNL calibration is however also affected by the EN and it 
is normally performed by blocking the input signal, thus 
replacing it by vacuum. It yields the variance 

- ^2fysNU' = 4aio(5^ic) + {^^lii) + (3) 

which is clearly a false estimation of the SNL due to the 
EN contribution. To reliably estimate the SNL, an inde- 
pendent measurement of the EN level must be executed 
(by blocking all input beams to the detectors). At very 
high LO powers, where the contribution from the EN is 
relatively low, it can be subsequently subtracted from the 
measurement outcome in ([3]) thus yielding 

VsNL = 4aio- (4) 

And the SNL for all powers can then be found by extrap- 
olation using the fact that the SNL scales linearly with 
the power of the LO. We note that such a calibration 
can be performed with any bright light beam provided 
that the homodyne detector is well balanced or the light 
beam is shot noise limited. There is no requirements for 
spatio-temporal matching with the signal beam (which is 
just the vacuum state). 

A traditional squeezing measurement is therefore car- 
ried out by recording the variance of the difference cur- 
rent (using a spectrum analyzer) and subsequently com- 
pare this to the shot noise level: 

^ ^ ((»! -i2Y)sX^ 
VsNL 

= (4"io(^^|) + + {5^l,)) (5) 

As clearly seen from this expression, the measurement 
of the degree of squeezing is hampered by the elec- 
tronic noise of the detectors; only if (Siln) + {Sili^) « 




FIG. 1: Schematic of the setup. The beam from the squeez- 
ing setup (OPA or Kerr based squeezing setup, see text) is 
interfering with the LO on a 50/50 beamsplitter BS, with 
a relative phase of 0. Dl and D2 are the analogue pho- 
todetectors. The photocurrents of the detectors ii and i2 are 
analyzed by the acquisition system. 

Instead of constructing the difference current, it is in- 
teresting to consider the covariance of the two detector 
outcomes: 

cov(zi,i2) = (ii«2) - («i)(«2) = (a|ai4a2) - 

(4ai)(4a2) = aloliSXlo) - (<5X|)], (6) 

Here it is assumed that the EN of the two detectors are 
not correlated i.e. {SieiiSien) — {Sieii){Siei2)- As seen 
from ([6]) the covariance coefficient is completely indepen- 
dent of the EN due to the time averaging of the data and 
the statistical independence of the noise of two photode- 
tectors. For the sake of clarity, we assume that the beam 
of the LO is shot noise limited i.e. (SX^^q) — (SX^^^). 
The analysis of the formula ^ suggests, that the sign 
of the covariance coefficient is determined by the statis- 
tics of the incoming light. Indeed, if the input state is 
squeezed i.e. {SX^) < (SX^^^), the covariance is pos- 
itive. In contrast, if the input state is coherent i.e. 
{SX'^) = (SX^^^) the covariance is equal to zero, and 
finally, if the input state exhibits classical excess noise 
i.e. (SX^) > (i5X^£,c), the covariance is negative. There- 
fore, the sign of the covariance coefficient allows one to 
carry out a qualitative characterization of the statistics 
of the measured quadrature. 

The covariance measurements can be also used to es- 
timate the degree of squeezing. However, in calibrating 
the shot noise level one must resort to the standard HD 
approach using the methods described above. Once this 
calibration is obtained for a specific pair of detectors it 
can be used for any subsequent squeezing measurements 
exploiting the covariance method. By comparing eq. Q 
for the shot noise with eq. ^ for the covariance we find 
the following expression for the degree of squeezing 

SH«|> = i-4£2Zi!i:M (7) 

VSNL 

This expression is seen to be independent on the elec- 
tronic noise of the detectors, and thus the correlation 
method yields an electronic noise free estimation of 
squeezing. The estimation relies, however, on a noise-free 
calibration of the SNL, that is, Vsnl must be measured 
accurately as mentioned above. 

Let us mention that in the derivation of the covariance 
coefficient (6), we did not perform a normal ordering pro- 
cedure of the creation and annihilation operators. On 
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the other hand, the joint detection operator, referred as 
a coincidence operator, which is traditionally introduced 
in discrete variable quantum optics, is defined through 
the second normally-ordered intensity moment [^. In 
order to find the link between the coincidences operator 
and the covariance coefficient, we rewrite the normally 
ordered coincidence operator and apply the commutation 
relations for a and a^: 

(: iiZ2 :) - (ii)(i2) = («i«2) - (i2) 

=^ (: 1112 :) = cov{ii,i2) + {ii){i2)- (8) 

Thus, the covariance coefficient, revealed from the di- 
rect product of two photocurrents, corresponds to the 
normally ordered coincidence operator (: iii2 :) with sub- 
tracted background of the " accidental coincidences" , pro- 
portional to (ii)(i2)- 

It is also useful to see how the covariance coefficient 
scales with linear losses. Let us first consider a situa- 
tion where the LO and the input signal experience the 
same linear losses rj. From a simple beamsplitter model 
it follows that the LO will loose in mean intensity and a 
portion of the vacuum fluctuations will be coupled to the 
input beam. We obtain the following expression for the 
attenuated covariance: 

cov,(ii,Z2) - (1 - vfalomlo) ~ (SXl)]. (9) 

This result is confirmed by the general rule, that the 
second normally ordered intensity moment scales with 
the square of the transmission, i.e. (: ii«2 '■)ri = (1 ~ 
7y)^(: iii2 :) 3- We now consider the scenario where the 
LO power is decreased but without affecting the signal 
state. In this case the covariance coefficient scales lin- 
early with the transmission, I — 77: 

cov,(*i,*2) = {I - vhloiiSXlo) ~ (SXl)]. (10) 

These linear and quadratic dependencies with respect to 
the transmission will be experimentally verified in the 
following sections. 

III. EXPERIMENT 

The idea described in the previous section was experi- 
mentally investigated in ref. [15j] for a very bright polar- 
ization squeezed state. In that realization, the LO was in 
the same spatial mode with the squeezed vacuum. There- 
fore, it was not possible to attenuate the two modes (the 
squeezed and the LO) independently. In the present pa- 
per we extend the previous analysis to include also the 
investigation of the squeezed vacuum prepared in a mode 
spatially separated from the LO, thus allowing indepen- 
dent control of the two modes. For completeness, we 
present in this and the following section both the results 
of the new experiment as well as the results of ref. [l^ . 

The setup for the bright squeezing experiment, based 
on the Kerr nonlinearity in optical fibers, is shown in 
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FIG. 2: The scheme of the polarization squeezing setup. The 
beam from the Cr''"'":YAG laser is split on the polarization 
beam splitter (PBS) after polarization rotation, performed 
by a half- wave plate (HWP). The birefringence compensator 
is used to provide a temporal overlap between two orthogo- 
nally polarized pulses after propagation through a single mode 
polarization maintaining fiber. The interference between the 
pulses is observed by placing a HWP and a PBS in front of two 
analogue photodetectors. The photocurrents of two detectors 
are further analyzed by the high-speed digitizer. 
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FIG. 3: The scheme of the OPA based setup. A frequency 
doubled Nd:YAG laser (fundamental wavelength at 1064 nm) 
is used to pump the bow type shaped OPA with a 10mm long 
PPTKP crystal. The squeezed vacuum field, generated by 
the OPA, is mixed with the fundamental wave of the pump 
laser on a 50/50 beam splitter (BS), used as a LO, and further 
detected by two analogue photodetectors (Detl, Det2). The 
relative phase between the squeezed vacuum field and the LO 
is controlled by a mirror on a piezo-stage (not shown). The 
photocurrents of two detectors are further analyzed by the 
high-speed digitizer. 



Fig. El [H, [H. A Cr^+iYAG laser generates pulses with 
a duration of 140 femtoseconds, at a repetition rate of 
163 MHz and with a central wavelength of 1497 nm. The 
laser pulse is first split between two orthogonal polar- 
ization modes as well as two different temporal modes 
using a "birefringence compensator" (see Fig. 2) and 
secondly launched into a 13.3 meter long single mode po- 
larization maintaining fiber. While propagating through 
the fiber, the orthogonally polarized pulses experience 
the Kerr effect, resulting in the squeezing of a certain 
quadrature in both polarization modes. The temporal 
overlap between two output pulses is provided by the 
" birefringence compensator" in front of the fiber as well 
as the locking loop based on a Stokes parameter mea- 
surement of a small portion (< 0.1%) of the fibre out- 
put. This resulted in a circularly polarized beam at the 

fiber output, mathematically described by ( 6*3 ) 7^ and 
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Si J — \^S2j — (where 5*1, 5*2 and S3 are the three 

Stokes parameters). The conjugate polarization opera- 
tors, which can exhibit polarization squeezing, are then 
found in the so called dark plane given by 5*1 — 52. The 
polarization squeezing is observed by a Stokes parame- 
ter measurement setup with the help of a HWP placed 
in front of a PBS. Then two analogue high-efficient pho- 
todetectors with quantum effiencies > 98% are placed 
in the output ports of the PBS and their photocurrents 
are analyzed by a high-speed digitizer. Such polarization 
squeezing is equivalent to a " vacuum" squeezing in a cer- 
tain polarization mode (here the left circular polarization 
mode) whereas the orthogonal polarization mode (right 
circular polarization mode) serves as the intrinsic local 
oscillator The relative phase between the LO and 

the squeezed "vacuum" is varied by rotating the HWP 
placed in front of the PBS in the Stokes measurement 
setup. Note that since the two polarization modes share 
the same spatial mode, independent control or attenua- 
tion of them is not straightforward. This limitation lead 
us to also use an alternative source of the squeezed light 
based on the optical parametric amplificr(OPA) which is 
described below. 

The OPA setup is shown in Fig. [31 We use a bow 
tie shaped cavity with a nonlinear crystal, pumped by 
the second harmonic (532 nm) of a continuous wave sin- 
gle mode field from a Nd:YAG laser. The OPA has a 
cavity round trip path length of 275mm and the trans- 
mittance of the output coupler (OC) is 10%. A 10mm 
long periodically poled KTP (PPKTP) crystal is used as 
a nonlinear medium inside the cavity, and the cavity is 
phase locked using the Pound-Drever-Hall locking tech- 
niques on counter propagating light field. The squeezed 
vacuum field at 1064 nm is measured by combining it on 
50/50 beam splitter with a LO drawn directly from the 
laser. The relative phase between the squeezed vacuum 
and the LO is controlled by a mirror on a piezo stage. 

In both experimental arrangements we recorded the 
raw data from the two detectors by digitizing an AC 
component of the photocurrents. For this purpose we 
used a high-speed digitizer with sampling rates of 2 • 
10®samples/s and 20-lO^samples/s and input bandwidths 
of 150KHz and 3MHz for the OPA and the Kerr based 
squeezing setup, respectively. Finally, a simple Mat- 
lab script evaluated the variance of the difference cur- 
rent (corresponding to standard HD) and the covari- 
ance between the two detector outcomes. This enabled a 
straightforward comparison between the two methods. 

By changing the phase (j> of the LO (using a half wave 
plate in Fig. [5] and mirror on a piezo stage in Fig. [3]), 
we were able to choose the measured quadrature 6X^ to 
be either squeezed or anti-squeezed. The measurement 
of the coherent state was taken with a shot noise limited 
beam, whilst the squeezed beam was blocked. 

The influence of the EN on the results of the mea- 
surements was studied by attenuating the intensity of 
the optical signal using neutral density (ND) filters. In 
case of strong attenuation, the produced photocurrent 
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FIG. 4: The covariance versus the phase of the LO at two dif- 
ferent powers of the LO (open and filled circles) for the OPA 
based squeezing setup. Positive covariance values correspond 
to the measurement of the squeezed quadrature, whilst the 
negative covariance values correspond to the measurement of 
the antisqueezed quadrature. 
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FIG. 5: The covariance at the maximally squeezed (open cir- 
cles) and the maximally antisqueezed quadrature (filled cir- 
cles) versus the transmittance of the LO (l-rj). 



will carry a relatively large contribution from the EN, 
whereas for no attenuation the EN plays only a minor 
role. Therefore through attenuation, the transition be- 
tween EN dominance and the optical noise dominance 
can be smoothly accessed. For the Kerr squeezed state 
both the LO as well as the squeezed state was attenuated 
whereas for the OPA experiment only the LO underwent 
the attenuation. 



IV. RESULTS AND DISCUSSION 

First we investigated the covariance as a function of 
the phase (j) of the LO. In Fig. 2] we present the results 
of such measurements using the OPA setup with two dif- 
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FIG. 6; Effect of the transmittance of the LO (l-J?) on the 
covariance at (a) the maximaUy squeezed and (b) the maxi- 
mally antisqueezed quadrature. In both plots we compare the 
covariance of a squeezed beam (red) to a calibration measure- 
ment using a coherent beam (black). 



ferent powers of the LO. Positive covariance corresponds 
to the measurement of squeezed quadratures, whilst the 
negative covariance corresponds to the measurement of 
the anti-squeezed quadratures. We see that the phase 
interval for which squeezing is witnessed is independent 
on the power of the local oscillator. With other words, 
despite the presence of a relatively large contribution of 
electronic noise in the detectors while using a weak LO, 
the presence of squeezing in the signal is still evident. 
A similar trend of the covariance was observed for the 
Kerr based squeezing source fl^. However, in that case 
the covariance modulus of the anti-squeezing quadrature 
were much larger due to the larger impurity of the Kerr 
squeezed state fis']. 

Correlation plots as the ones in Fig. [3] were made for 
various different powers of the LO, and for each real- 
ization we recorded the maximum correlation and anti- 
correlation. The results are shown in Fig.[5]where the co- 
variance is plotted against the power of the LO. It clearly 
demonstrates the linear dependence as predicted by equa- 
tion (jlOp . We performed similar measurements for the 
Kerr squeezed state. This time, however, the squeezed 
mode as well as the LO mode were attenuated simul- 
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FIG. 7: The degree of squeezing (triangles) and antisqueez- 
ing (circles) revealed from the correlation analysis versus the 
transmittance of the LO (I-77) in the OPA based squeezing 
setup. The relatively large error bars are mostly due to sta- 
tistical fluctuations. 

taneously, which according to eq. ([9]) yields a quadratic 
dependence in l-rj. This was experimentally confirmed 
by the measurements as summarized in Fig. [6l 

In the next step of the experiment we wanted to per- 
form an electronic noise free measurement of the degree 
of squeezing. For this purpose we first calibrated the 
shot noise level. It was important that this calibration 
was not affected by the electronic noise of the detectors 
as this otherwise would yield a false estimation of the 
SNL and a corresponding overestimation of the degree 
of squeezing. The calibration was therefore established 
at very high powers of the LO, making the first term in 
eq. ([3|) very large and the following EN terms negligible. 
Using the fact that the shot noise scales linearly with 
power, we could easily determine the shot noise at all 
powers from the high power measurement data. These 
calibrations (one for the OPA experiment and one for 
the Kerr experiment) were then subsequently used for 
establishing the degree of squeezing. 

From the correlation measurements, presented in Fig. [5] 
and[6l and the corresponding shot noise calibrations, we 
use eq. ([7]) to compute the degree of squeezing. The re- 
sults are summarized in Fig.[7]for the OPA squeezer and 
Fig. [H] for the Kerr squeezer. We clearly see from Fig. [7| 
that the measured degree of squeezing is almost constant 
as the LO power is decreased although the relative con- 
tribution from the EN becomes large. This should be 
compared with traditional homodyne detection of the de- 
gree of squeezing, the measurements of which are shown 
in Fig. [HI Here we see the expected trend for homodyne 
detection, that the measured degree of squeezing goes 
down as the relative EN increases. This is indeed the 
main conclusion of this paper; by using the correlation 
method, it is possible to reliably estimate the degree of 
squeezing independent of the EN (even though the LO 
is very dim) in contrast to the standard homodyne mea- 
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FIG. 8: Measurement along the squeezed quadrature in the 
OPA based squeezing setup using homodyne detection. The 
plot shows, how the degree of squeezing depends on the trans- 
mission of the LO (l-ri). 
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surement approach which is affected by electronic noise. 
This also holds true for a LO free setup where a very dim 
squeezed beam is directly measured (using a beam split- 
ter and two detectors). In this case, however, the signs 
of the correlations will swap and the squeezing will be 
computable from the expression S = ^'^°v^l^^^^ ~ 1 ^M- 

Similar investigations were carried out for the Kerr 
squeezer as described in ref.(l5l| and summarized in 
Fig. [HI in Fig. [9l\. and B we plot the measured degree 
of squeezing versus the transmission using homodyne de- 
tection and correlation measurements, respectively. The 
insets show magnified parts of the plots. Similar con- 
clusions can be drawn from these plots as from Figs. [7] 
and [51 The correlation method is capable of measuring 
squeezing at lower LO powers than for the homodyne 
detectors. In these experiments, however, the attenua- 
tion was applied collectively on both modes which means 
that the actual degree of squeezing was degraded as the 
transmittivity goes down. 

We should emphasize that the "removal" of electronic 
noise is a result of the averaging procedure after detec- 
tion as mentioned in section II. Every single measure- 
ment data is largely affected by electronic noise which 
means that the correlation method is not yielding elec- 
tronic noise free data " in real time" but only after several 
measurements of an ensemble and subsequently averag- 
ing. Moreover, the single measurement does not give 
information about the amplitude of the signal and thus 
cannot in general be used in feedforward based quantum 
information protocols. 



CONCLUSIONS 
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FIG. 9: Experimental dependence of the amount of squeezing 
on the transmission (a) for homodyne detection and (b) 

for the covariance method. The solid lines represent the the- 
oretical prediction without accounting for the EN. The hori- 
zontal dashed line at zero marks the border between squeezing 
and classical excess noise. 



In conclusion, we have studied the implementation of 
a correlation measurement technique on two different 
sources of squeezed light generated from an optical para- 
metric oscillator and from an optical fiber. We showed 
that this method can be used to estimate the degree of 
squeezing independent of the amount of electronic noise 
in the detectors provided that their noises are uncorre- 
lated. This is in stark contrast to using the standard 
homodyne detection to estimate the degree of squeezing, 
where electronic noise might affect the measurements, es- 
pecially when the LO power is low. 

The correlation measurement strategy is particular 
useful for measuring the degree of squeezing of a squeezed 
light that is neither matchable to any local oscillator 
modes nor has a high coherent excitation. For such cases, 
the electronic noise becomes a critical issue which can be 
overcome by the correlation measurement method. 
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